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ABSTRACT 

We present J-band spectroscopy of passive galaxies focusing on the Na 1 doublet at 1.14 fim. 
Like the Na 10.82 pm doublet, this feature is strong in low-mass stars and hence may provide a 
useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, 
we find that Nal 1.14 pm increases steeply with increasing velocity dispersion, cr, and for 
the most massive galaxies (cr > 300km s“^) is much stronger than predicted from synthetic 
spectra with Milky-Way-like IMFs and solar abundances. Reproducing Nal 1.14pm at high¬ 
er likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of 
both. Using the NaD line to break the degeneracy between IMF and abundance, we infer 
[Na/H] Ri h- 0.5 and a steep IMF (single-slope-equivalent a; « 3.2, where x = 2.35 for Salpeter), 
for the high-cr galaxies. At lower mass (cr = 50-100km s“^), the line strengths are compatible 
with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where 
strong gravitational lensing masses favour MW-like IMFs. Like the high-cr sample on average, 
these galaxies have strong Nal 1.14 pm; taken in isolation their sodium indices imply bottom- 
heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum- 
fitting approach, applied to the high-cr sample, recovers an IMF less heavy than Salpeter, but 
under-predicts the Nal 1.14pm line at the 5a level. We conclude that current models struggle 
to reproduce this feature in the most massive galaxies without breaking other constraints, and 
caution against over-reliance on the sodium lines in spectroscopic IMF studies. 

Key words: galaxies: stellar content — galaxies: elliptical and lenticular, cD 


1 INTRODUCTION 

Because the luminosity-vs-mass relation for stars is steeper than the 
distribution function of their masses, dwarf stars contribute very lit¬ 
tle to the integrated light from an old stellar population. However it 
has long been realised that some gravity-sensitive spectral features, 
which become extremely strong in late M dwarfs, provide a route 
to constraining the low-mass stellar content, and hence the initial 
mass function (IMF), in unresolved galaxies (e.g. Spinrad & Tay¬ 
lor 1971; Whitford 1977; Carter, Visvanathan & Pickles 1986). The 
Na I doublet absorption at 0.82 pm is among the strongest and most 
readily accessible of these dwarf-star indicators. 

Recent work has advanced this method using measurements 
with red-sensitive CCDs, large spectroscopic surveys, and im¬ 
proved near-infra-red (IR) instrumentation (e.g. van Dokkum & 
Conroy 2010; Spiniello et al. 2012; Conroy & van Dokkum 2012b; 
Smith, Lucey & Carter 2012; La Barbera et al. 2013). In parallel, 
more extensive stellar libraries (e.g. Rayner et al. 2008) and so¬ 
phisticated spectral synthesis models (e.g. Conroy & van Dokkum 
2012a, CvD12a; Vazdekis et al. 2012) have led to better quantifica¬ 
tion of the IMF sensitivity of various features, and the degeneracies 
with element abundance variations. 

These studies found that the Na I 0.82 pm doublet observed 


in massive early-type galaxies is much stronger than would be ex¬ 
pected from an old stellar population with an IMF like that in the 
Milky Way, and with solar element abundances (van Dokkum & 
Conroy 2010; Spiniello et al. 2012; Conroy & van Dokkum 2012b; 
Fererras et al. 2013). The observations can be reproduced by adopt¬ 
ing either a bottom-heavy IMF or a highly-enhanced sodium abun¬ 
dance, or a combination of both. As an example, for NGC 4621, the 
full-spectrum fitting analysis of Conroy & van Dokkum (2012b) 
requires an IMF with a mass-to-light ratio larger by a factor 
Q= 1.95±0.05, compared to the Milky Way IMF, and a sodium 
abundance [Na/H] =-l-0.84±0.04 (their figure 3). However, the re¬ 
sults are sensitive to the spectral regions used in the fit: if Na I 
0.82 pm is excluded, they derive a fa 1.2 for this galaxy, not much 
in excess of a Kroupa (2001) IMF (a = 1), while if the Ca II triplet 
is also excluded, a heavy IMF is again recovered (their figure 12). 

Since Nal 0.82pm plays an important role in the evidence 
for non-standard IMFs in massive galaxies, and because alter¬ 
native model sets disagree in predicting the Nal 0.82pm and 
NaD line behaviour (Spiniello, Trager & Koopmans 2015), it is 
worthwhile not only to explore other indicators based on different 
atomic/molecular species (e.g. Spiniello et al. 2014), but also to test 
models against a wider set of observed sodium features. 

In this Letter, we focus on the J-band Nal doublet absorp- 
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tion at 1.14/im. This line arises from atoms in the same energy 
level as Nat 0.82/rm, and hence should have similar “fundamen¬ 
tal” temperature and pressure dependence, although blending with 
other features inevitably affects the behaviours differently in prac¬ 
tice (e.g., Nat 0.82/rm coincides with a strong TiO band). As ex¬ 
pected, Nat 1.14 /rm is very strong in cool dwarf stars, and weak in 
giants (Rayner et al. 2008), yielding good predicted sensitivity to 
the IMF (Conroy & van Dokkum 2012a). No published work has 
yet analysed Nat 1.14/rm in early-type galaxies, most IR studies to 
date having been made at longer wavelengths (Silva, Kuntschner & 
Lyubenova 2008; Marmol-Queralto et al. 2009; Cesetti et al. 2009). 
For very nearby galaxies, accurate measurement of Nat 1.14 ^m is 
hampered by telluric water absorption at 1.11-1.16/rm. Flere, we 
study two samples of galaxies at modest distance {z = 0.02-0.06), 
so that Nat 1.14/rm is redshifted into a cleaner atmospheric win¬ 
dow. We use rest-frame stacking to build composite spectra, further 
suppressing the sky and telluric noise contributions. 

As this is the first analysis of the Nat 1.14 ^m feature in this 
context, our treatment is exploratory in nature, and intended to cap¬ 
ture the “first order” behaviour and to distinguish IMF from sodium 
abundance effects, rather than to provide an exhaustive treatment of 
all possible parameter degeneracies. 

Our samples, observations and data treatment are described in 
Sectionl^and the results are presented in Section|^ We summarize 
our conclusions in Section|4] 


2 DATA 

We combine J-band data from two observational programmes to 
measure Na 11.14 /rm over a wide range in galaxy mass, and derive 
NaD from optical survey spectroscopy for the same galaxies. 

2.1 Subaru Coma sample 

We obtained J-band spectra for 84 red-sequence galaxies in the 
Coma cluster, with velocity dispersions a = 50-300 km s“^, using 
FMOS, the IR Fiber Multi-Object Spectrograph at the 8.2m Subaru 
Telescope (Kimura et al. 2010). The observations were made on the 
nights of 2014 April 9-10, in the “J-long” configuration, with spec¬ 
tral resolution R « 2400. The target galaxies, fibre configurations 
and observing strategy were the same as those used in the Y-band 
(“J-short”) observations reported by Smith et al. (2012b). Four fi¬ 
bre configurations were observed in cross-beam-switching mode, 
with fainter galaxies repeated on more configurations, for a total 
integration of up to 10.5 h per target. 

Initial data reduction was performed using the standard FMOS 
FIBRE-PAC pipeline (Iwamuro et al. 2012), yielding extracted, cal¬ 
ibrated spectra for each fibre. A correction for telluric absorption 
was made using stars observed in the field simultaneously with 
the galaxy targets. Because the FMOS OH suppression mask re¬ 
moves significant parts of each spectrum, rest-frame stacking is 
essential to derive clean and continuous spectra from these obser¬ 
vations. We constructed two independent composite spectra from 
the FMOS data, according to literature velocity dispersion mea¬ 
surements compiled by Smith et al. (2012a). The low-cr stack in¬ 
cludes 32 galaxies with literature velocity dispersions in the range 
50 < CT < 100kms“^, while the higher-cr stack has 52 galaxies with 
100 < a < 300kms“^. By bootstrap resampling the input spectra 
used in the stacking process, we computed a set of 100 alternative 
realisations of the composite spectra, which are used to estimate 
the errors, including systematic variation among galaxies within 


the stack. The signal-to-noise ratio (S/N) estimated from the boot¬ 
straps is 100 (lower a) and 200 (higher cr) at 1.14 ^m. 


2.2 VLT SNELLS lens-candidate sample 

Very massive galaxies (a > 300 km s“^) are not well represented 
in the Coma sample. Since these objects are of particular inter¬ 
est, given the evidence for non-standard IMFs in the most massive 
galaxies, we use a second dataset to sample this regime. 

The observations for the high-mass sample were obtained for 
the “SNELLS” (SINFONI Nearby Elliptical Lens Locator Survey) 
project, a search for low-redshift strong-lensing galaxies among the 
massive early-type population (Smith, Lucey & Conroy 2015). The 
targets were drawn from the Six-degree Field (6dF) Peculiar Veloc¬ 
ity Survey (Campbell et al. 2014), and the Sloan Digital Sky Sur¬ 
vey (SDSS, Data Release 7, Abazajian et al. 2009), with primary 
selection on velocity dispersion, with a > 300kms“^ and redshift 
2 < 0.055. Galaxies with strong emission lines were excluded. We 
observed 26 candidate lenses in the J-band with the SINFONI 
integral-field spectrograph at the ESO 8.2 m VLT UT4. The largest 
image scale was used, to provide a field-of-view 8x8 arcsec^. The 
spectral resolution is i? « 2000. The observations were optimized 
for efficiency in finding background line emitters, with two 600 sec 
exposures obtained per galaxy, offset spatially by 2.3 arcsec. 

The spectrum for each galaxy was extracted from all spatial 
pixels with (absolute) total flux >50 per cent of the maximum, af¬ 
ter subtracting the two datacubes to remove the background. This 
corresponds to an aperture of ~l-2 arcsec^ in most cases. A telluric 
absorption correction was determined by fitting water and molec¬ 
ular oxygen transmission functions directly to the extracted spec¬ 
trum, after dividing out a representative passive galaxy model spec¬ 
trum. As with the EMOS data, we combine the SNELLS spectra in 
the rest frame to derive a clean composite spectrum, and compute 
bootstrap realisations to quantify the errors. The derived S/N is 
~350 A“^ at 1.14 ^m. The effective velocity dispersion, measured 
from the final J-band stack, is 297±7 km s“^, slightly smaller than 
expected given the SDSS/6dF measurements. 

Among the SNELLS galaxies are two which were confirmed 
as strong lensing systems (SNL-1 and SNL-2) by Smith et al. 
(2015). The lensing analysis shows that these two galaxies have 
mass-to-light rations similar to what is expected for a MW-like 
IMF Specifically, the mass-excess factors, relative to the Kroupa 
IMF, are a=L20±0.13 and a = 0.94±0.17; hence bottom-heavy 
models are strongly disfavoured in these cases. 


2.3 Optical spectra 

To aid in interpreting the Nat L14/im observations, we also use 
optical spectra to measure the NaD line at 0.58 ^m. The data are 
drawn from the 6dF survey (for most of the SNELLS galaxies) and 
SDSS (for the Coma sample, and the remaining SNELLS galaxies). 
We created composite spectra and bootstrap realisations from the 
same galaxies as used in the FMOS and SNELLS stacks. 

We note that the aperture sizes probed by the various datasets 
are not identical: the fibre sizes are 1.2 arcsec, 3 arcsec and 
6.7 arcsec in diameter, for FMOS, SDSS and 6dF respectively, 
while the effective aperture for the SNELLS spectra is 1-2 arcsec. 
This is especially relevant when comparing the NaD measure¬ 
ments, taken from large apertures (SDSS and 6dF), against the J- 
band spectra obtained in small apertures from SNELLS and FMOS. 
Although the SNELLS spectra are derived from IFU data, our 
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Figure 1. Upper panel: composite J-band spectra for early-type galaxies in three bins of velocity dispersion. Lower left: stacked SDSS and 6dF spectra for 
the NaD region, using the same galaxies. Lower centre: zoom on the Nal 1.14/rm doublet region in the J-band. Lower right: zoom on the weak Nal line at 
1.27 ^m. In each panel, the thick black lines show the stacked galaxy spectra, while the orange crosses are the individual flux points for two confirmed lenses 
SNL-1 and SNL-2 from Smith et al. (2015). Thin lines show the predictions from illustrative models from Conroy & van Dokkum (2012a). The observed and 
model spectra have been smoothed to match the velocity dispersion of the SNELLS sample. 


observing strategy is not suited to extracting spectra at larger aper¬ 
ture to match the SDSS/6dF spectra. However, we can estimate the 
likely effect of aperture mismatch using NaD measured in a range 
of aperture sizes in data cubes from the CALIFA survey (Sanchez 
et al. 2012). For 22 elliptical galaxies in the CALIFA second data 
release, the NaD line strength is reduced by 9 per cent on average, 
with a scatter of 5 per cent, when the aperture size is increased by a 
factor of five. Correcting NaD upwards, to account for the smaller 
apertures in SNELLS and FMOS, does not affect our conclusions, 
as we show in Section|L2] 


3 RESULTS 

We present the results both through qualitative comparison of ob¬ 
served spectra against synthetic spectra, and through a more quanti¬ 
tative analysis using Lick-style line indices. In both approaches, we 
use improved versions of the CvD models, as described by Conroy, 
Graves & van Dokkum (2014). 

3.1 Spectral comparison 

The stacked spectra from the FMOS and SNELLS samples, to¬ 
gether with the corresponding NaD stacks, are shown in Eigure[T] 
with illustrative model spectra for comparison. The EMOS stacks 
and the models have been smoothed to match the velocity disper¬ 
sion of the SNELLS galaxies. Prominent spectral features in the 
J-band include the Nal L14^m doublet itself and the hydrogen 
Pa P line at 1.28 /rm, as well as many other heavily blended metal 
lines (see Rayner et al. 2008 for identifications). 


As a fiducial comparison, we consider the 13.5 Gyr 
model with Kroupa IMF and abundance pattern [a/H] = -F0.3; 
[Fe/H] = [Na/H] = 0. This model clearly under-predicts both the 
Nat L14/im and the NaD absorption in the SNELLS sample. 
The discrepancy is smaller in the FMOS stacks. The predicted 
Nat 1.14/rm strength can be increased by adopting either a more 
bottom-heavy IMF or larger sodium abundances, but NaD is lit¬ 
tle affected by the IMF according to the CvD models. Hence some 
enhancement in [Na/H] is essential to match the SNELLS sample. 
This is supported by the noticeable presence of a weak Na I line at 
1.27 fim, which only becomes visible in models with high [Na/H]. 
Taken individually, the two gravitational lenses, SNL-1 and SNL-2, 
both have spectra consistent with the SNELLS high-cr composite, 
and in particular follow the behaviour of the stack in the NaD and 
Nat 1.14/rm regions. 

Eigurel^shows the ratio between the observed SNELLS stack 
and the fiducial model, compared to ‘response spectra’ show¬ 
ing the effect of some relevant stellar population parameters. The 
observed-vs-fiducial ratio spectrum confirms the excess absorption 
in the sodium lines, relative to the model. Only the [Na/H] and IMF 
responses have a strong impact on the Nat L14^m absorption. 
While either sodium enhancement or a bottom-heavy IMF could 
reproduce Nal L14/rm, changes to these parameters lead to dif¬ 
ferent predictions elsewhere, most notably at NaD, but also in the 
weak sodium lines at 1.23 fim and 1.27 /im, and the potassium dou¬ 
blets at 1.17 fim and 1.25 fim. Visual inspection of the Na I 1.14 ^.m 
and NaD effects suggests a combination of sodium enhancement 
and bottom-heavy IMF is required for the high-a galaxies, if the 
combined effect of other parameters (age, etc) can be neglected. 
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Figure 2. IMF and abundance effects on the J-band spectra and the Na D region. The black line shows the ratio of the observed SNELLS stack versus a fiducial 
model (with [cs/H] =+0.3, [Fe/H] = [Na/H] = 0, 13.5 Gyr age and Kroupa IMF). The blue and red lines show the response spectra for sodium enhancement and 
IMF variation respectively. The effects of some other relevant parameters are shown below. The age response is shown both for the Kroupa IMF (solid) and 
for a; = 3 (dashed). Because NaD is intrinsically much stronger than the IR features, the ratio spectra in the NaD region have been scaled by a factor of 0.07 
relative to the J-band for clearer presentation. Shaded regions indicate the feature and pseudo-continuum definitions of the Na 11.14 fim and Na D indices used 
in Section[T2] 


3.2 Line index comparison 

To summarize the behaviour of the strong sodium lines, we have 
measured Nat 1.14^m and NaD using absorption indices. For 
NaD, we use the Lick definition from Trager et al. (1998), while for 
Nal 1.14^m, we adopt a feature band of 1.1367-1.1420/rm, with 
pseudo-continua at 1.1337-1.1367 ^m and 1.1420-1.1450 ^m (see 
Figure O. All of the results are robust to small changes to these 
bandpasses. ForFMOS, the indices were measured after smoothing 
to match the velocity dispersion of the SNELLS stack. The index 
errors were derived from the bootstrap stacks, and hence explicitly 
include the galaxy-to-galaxy scatter within each sample. 

Figure [3 shows the measured indices, compared to a model 
grid spanning a range of [Na/H] and different IMFs. The in¬ 
dices were measured on the model spectra after smoothing to 
(j = 300kms“^ to match the data. Note that the model grid is for 
[q:/H]=+ 0.3, which is appropriate for the SNELLS galaxies, but 
somewhat high for the EMOS stacks. (From optical spectra, the 
FMOS samples have average [Mg/Fe] = 0.15 and 0.24.) 

The index-index diagram confirms and quantifies the results 
seen in the spectral comparisons: Nal 1.14/rm increases strongly 
with galaxy mass, deviating far from the predictions for solar 
[Na/H] and MW-like IMFs. NaD also increases strongly with cr, 
but not steeply enough to match a constant IMF model. Instead, 
the diagram implies a combination of both increasing [Na/H] and 
increasing dwarf-star content. For the SNELLS stack, the grid sug¬ 
gests an IME slope a:«3.1 [where the Salpeter (1955) slope is 
2.35] and [Na/H] « +0.5. Here, x should be interpreted as a single- 
slope-equivalent value, since the detailed shape of the low-mass 
IME is not constrained. Applying the CALIFA-derived aperture 
trend to correct NaD to the same aperture as used in the IR mea¬ 
surement, the inferred IMF slope would be reduced slightly, to 


X « 3.0, and the sodium abundance increased to [Na/H] « +0.6. In 
the low-cr FMOS stack, comparison with a model grid at lower 
[a/H] (not shown in the figure) leads to a derived [Na/H] «-0.3, 
and an IMF consistent with either Kroupa or Salpeter. The run of 
[Na/H] is thus similar to that found from stacked SDSS spectra by 
Conroy, Graves & van Dokkum (2014), who report [Na/H] =-0.24 
at (cr) =87kms“^ rising to [Na/H] =+0.43 at (cr) = 295 kms“’^. 

NaD is a resonance line, subject to contamination by inter¬ 
stellar absorption. Correcting for this would decrease the derived 
[Na/H] and increase the required x. However leong et al. (2013) 
studied “NaD excess” objects from SDSS and concluded that the 
interstellar contribution is negligible for ordinary early-type galax¬ 
ies. Note that for steep IMFs, assuming ages younger than the nom¬ 
inal 13.5 Gyr leads to weaker predicted absorption at Nat 1.14^m 
(and at NaD), and hence larger derived x and [Na/H]. 

Figure [3 also shows the sodium indices for all of the individ¬ 
ual SNELLS galaxies, highlighting the lenses SNL-1 and SNL-2. 
There is a wide spread in Nat 1.14/cm among the sample, span¬ 
ning the full range of model IMFs, at [Na/H] =0.3-0.5. As ex¬ 
pected from the similarity of their spectra to the SNELLS compos¬ 
ite, SNL-1 and SNL-2 have line strengths consistent with the prop¬ 
erties of the stack, yielding [Na/H] «+0.4 and a: = 3.0-3.5. Thus 
despite having mass-to-light ratios consistent with a MW-like IMF, 
these galaxies have sodium indices characteristic of very steep IMF 
slopes, at least if interpreted using otherwise “vanilla” versions of 
the CvD models (old, a-enhanced, other abundances solar). 


4 SUMMARY AND DISCUSSION 

The immediate result of this paper is that the Nat 1.14/am doublet 
is very strong in the most massive early-type galaxies. Synthetic 
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models with Milky-Way-like IMF and solar [Na/H] fail to match 
the observed absorption. This behaviour is similar to that seen in 
the optical Nal 0.82 /rm doublet, and suggests that the optical re¬ 
sults are not due to the specific spectral “environment”, e.g. contri¬ 
butions from the coincident Ti O band. Comparison with the CvD 
models suggests that reproducing the observed Nal 1.14^m re¬ 
quires high sodium abundances or bottom-heavy IMFs in massive 
galaxies. We show that consistent results for Na I 1.14 /rm and Na D 
can be obtained using a combination of both effects, with other pa¬ 
rameters of the model held at plausible values. For galaxies with 
a > 300 km s“^, we derive [Na/H] Ri -1-0.5 and an IMF with (single- 
slope-equivalent) a;«3. Two galaxies in our sample have strong- 
lensing masses compatible with an Milky-Way-like IMF, but their 
sodium indices imply steep IMFs as in the high-a stack. If the IMF 
really was a single power-law, the spectroscopic M/L exceeds the 
leasing constraint by a factor of q ~ 4. Alternative IMF shapes, 
can reduce but not eliminate the discrepancy (e.g. fig. 5 of Smith 
2014). For example the broken power-law assumed by La Barbera 
et al. (2013) would yield a 2.5 

We stress that the mass discrepancy applies specifically to re¬ 
sults from the sodium indices, and not the general approach of fit¬ 
ting CvD models to constrain the IMF. For example, we have ap¬ 
plied the full spectral fitting machinery of Conroy et al. (2014) to 
an extended SNELLS composite spectrum, including wider optical 
coverage from 6dF/SDSS as well as SINFONI data in both J and H. 
The best-fitting model has a modest mass excess factor of q « 1.4, 
and [Na/H]«-1-0.7, and fits well over most of the spectral range. 
This model has much more freedom than the simple grid shown 
in Figure [3 with additional parameters describing the abundance 
pattern, two variable low-mass IMF slopes, and various nuisance 
effects. The greater flexibility allows the fit to reach closer to the 
observed line strengths than would be expected from the grid, given 
the recovered Salpeter-like IMF, but the model still under-predicts 
the measured Nal 1.14 ^m index by ~5 (t. 

Our results suggest that the CvD models do not yet reproduce 
the strong Nal 1.14 fim in the most massive galaxies without vio¬ 
lating other constraints (lensing masses, full-spectral fit). Such dis¬ 
crepancies in specific spectral features underline the advantage of 
combining many gravity-sensitive indicators (e.g. Na I, Ca II, Fe H) 
in deriving spectroscopic constraints on the IMF. They may also 
help to identify limitations in the current models, and the corre¬ 
sponding opportunities for future improvements to the method. 
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